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£4 ' Abstract 

bp: 

Coalescing compact star binaries are expected to be among the strongest sources 
of gravitational radiation to be seen by laser interferometers. We present calculations 
of the final phase of inspiral of equal mass irrotational neutron star binaries and 
strange quark star binaries. Six types of equations of state at zero temperature 
are used - three realistic nuclear equations of state of various softness and three 
different MIT bag models of strange quark matter. We study the precoalescing 
stage within the Isenberg- Wilson-Mathews approximation of general relativity using 
a multidomain spectral method. The gravitational-radiation driven evolution of the 
binary system is approximated by a sequence of quasi-equilibrium configurations at 
fixed baryon number and decreasing separation. We find that the innermost stable 
circular orbit (ISCO) is given by an orbital instability for binary strange quark stars 
and by the mass-shedding limit for neutron star binaries. The gravitational wave 
frequency at the ISCO, which marks the end of the inspiral phase, is found to be 
~ 1100 — 1460 Hz for two 1.35 Mq irrotational strange stars described by the MIT 
bag model and between 800 Hz and 1230 Hz for neutron stars. 
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1 Introduction 



Coalescing neutron star binaries are considered among the strongest and most 
likely sources of gravitational waves to be seen by VIRGO/LIGO interferome- 
ters 11, 23, ij]. Due to the emission of gravitational radiation, binary compact 
stars decrease their orbital separation and finally merge. We separate the evo- 
lution of a binary system into three phases : point-like inspiral where orbital 
separation is much larger than the neutron star radius, hydro stationary inspi- 
ral where orbital separation is just a few times larger than the radius of the 
neutron star so that hydrodynamics play an important role, and merger in 
which the two stars coalesce dynamically. Gravitational waves emitted during 
last orbits of inspiral phase and the merger phase could yield important infor- 
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mation about the equation of state (EOS) of dense matter 1141. 13 _ 
Up to now, all relativistic calculations (except those of [29], J 25|) of the last 
orbits of inspiral phase have been done for the simplified EOS of dense matter - 
the polytropic EOS. In the paper we present results of our studies on late stage 
of inspiral of binary systems containing equal mass compact stars described 
by different realistic EOS of dense matter. The calculations are performed in 
the framework of Isenberg- Wilson-Mathews approximation to general relativ- 
ity (see Ref. [3j for a review). We consider binary systems consisting of two 
identical stars. We choose the gravitational mass of each star to be 1.35 M Q 
at infinite separation in order to be consistent with recent population syn- 
thesis calculations 0, [l9| and with the current set of well-measured neutron 
star masses in relativistic binary radio pulsars (2^, [h|. We assume that the 
velocity flow in the stellar interiors is irrotational (the fluid has zero vorticity 
in the inertial frame) since the viscosity of neutron star matter (or strange 
star matter) is far too low to ensure synchronization during the late stage of 
the inspiral 0, 24] . 



2 Equations of state and stellar models 



In Fig. 1 we show gravitational mass versus stellar radius for sequences of static 
compact stars. We limit ourselves to neutron stars consisting of nucleons and 
hyperons and strange quark stars described by the MIT bag model. Depending 
on the EOS we obtain the radius of a 1.35 M & star in the range 10-14 km. We 
perform calculations for three nuclear EOS of dense matter based on modern 
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Fig. 1. Gravitational mass M versus stellar radius R for sequences of static neutron 
stars described by three different nuclear equations of state (solid lines) and strange 
quark stars described by different MIT bag models (dashed lines). 

many-body calculations. The differences in the M-R relation for stellar models 
of neutron stars (solid lines) shown in Fig. 1 reflect the uncertainties in the 
existing theories of the interactions in nuclear matter. We consider one soft 
(BPAL12, U) and one stiff (AkmalPR, [2]) EOS of matter composed of 
nucleons, electrons and muons. We considered also one EOS in which hyperons 
are present at high densities (GlendNH3, jiol). The neutron star crust is 
described by means of a realistic EOS obtained in the many-body calculations 
(see jij for details). We note that the equations of state of nuclear matter can 
be considered static, as the timescales for bulk processes such as viscosity is 
much longer than the last stages of inspiral considered here. 



Strange stars are currently considered as a p ossible alternative to neutron 
stars as compact objects (see e.g. [33j, [271, ll8[ and references therein). Typi- 
cally, they are modeled with an EOS based on the MIT-bag model in which 
quark confinement is described by an energy term proportional to the vol- 
ume (e.g. [H, [2l|) in which quark confinement is described by an energy term 
proportional to the volume [15|. There are three physical quantities entering 
the MIT-bag model: the mass of the strange quarks, m s , the bag constant, B, 
and the strength of the QCD coupling constant a. In the framework of this 
model the quark matter is composed of massless u, d quarks, massive s quarks 
and electrons. Strange stars are self-bound objects, having high density (in 
the range ~ 3 — 6.4 [10 14 g/cm 3 ] 18]) at the surface. Dashed lines in Fig. 1 
correspond to sequences of static strange quark stars described by three dif- 
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Fig. 2. Binding energy as a function of gravitational wave frequency (twice the 
orbital frequency) along evolutionary sequences of irrotational binaries. The solid 
lines denote neutron stars, the dashed one strange quark stars and the dotted line 
point-mass binaries in the third order Post-Newtonian approximation 0] ■ Each evo- 
lutionary sequence finishes at the innermost stable circular orbit. The diamonds 
correspond to dynamical orbital instability while circles to the mass-shedding limit. 

ferent sets of parameters of the MIT-bag model: SQS56 - the standard MIT 
bag model: m s c 2 = 200 MeV, a = 0.2, B = 56 MeV/fm 3 ; SQSB60 - the 
simplified MIT bag model with m s = 0, a = 0; B = 60 MeV/fm 3 ; SQSB40 - 
the "extreme" MIT bag model (relatively low strange quark mass and B but 
high a) : m s c 2 = 100 MeV, a = 0.6, B = 40 MeV/fm 3 . 



3 Numerical method 



In the inspiral phase, since the timescale of orbital shrinking due to the emis- 
sion of gravitational waves is longer than the orbital period, one may consider 
a binary system to be in a quasi-equilibrium state (helical Killing vector ap- 
proximation). For each EOS, we construct so called an evolutionary sequence 
by calculating a sequence of quasi-equilibrium configurations with fixed baryon 
mass and decreasing orbital separation. We assume that the timescale in which 
the neutron star adjusts to the current gravitational potential is much shorter 
that the timescale of the orbital shrinking. The present computations of close 
binary neutron star or strange star systems rely on the assumption of irro- 
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tational flow of the fluid and a conformally flat spatial 3-metric (Isenberg- 
Wilson-Mathews approximation). In order to calculate the last orbits of inspi- 
ral phase of binary compact stars we use a numerical code which solves the five 
coupled, nonlinear, elliptic equations for the gravitational field, supplemented 
by an elliptic equation for the velocity potential of irrotational flows 



sec 
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for a discussion on different boundary conditions in the case of strange stars 
and neutron stars). The code has been already used successfully for calculating 
the final phase of the inspiral of compact stars [iol, 31, 32, 0, 25]. This code is 
built upon the C++ library LORENE (http : / /ww w . lorene . obspm . f r[ ). The 
complete description of the resulting general relativistic equations, the whole 
algorithm, as well numerous tests of the code can be found in 17]. Additional 



tests have been presented in Sect. Ill of [32 



4 Results 



In Fig. 2 we show the evolution of equal mass binary neutron stars (solid lines) 
and strange stars (dashed lines) having total gravitational mass 2.7M at in- 
finity. The binding energy Eu n d is defined as the difference between Madm 
(Arnowitt-Deser-Misner mass - the total mass-energy of a binary system [32j |) 
and the total mass of the system at infinite separation. This energy is equal to 
total energy emitted by a binary system in gravitational waves. The frequency 
of gravitational waves is twice the orbital frequency. Comparison of our nu- 
merical results with 3rd order PN point masses calculations 0] reveals a good 
agreement for small frequencies (large separations). The deviation from PN 
curves at higher frequencies (smaller separation) is due to hydrodynamical ef- 
fects, which are not taken into account in the PN approach. A turning point of 
-E-bind along an irrotational evolutionary sequence indicates the orbital (dynam- 
ical) instability fl6| . This instability originates both from relativistic effects 
and hydrodynamical effects. In the case where no turning point of -Ebind occurs 
along the sequence, the mass-shedding limit (Roche lobe overflow) marks the 
end of the inspiral phase of the binary system, since recent dynamical calcula- 
tions for 7 = 2 polytrope have shown that the time to coalescence was shorter 
than one orbital period for configurations at the mass-shedding limit (i.e. see 
(j^l ). Thus the physical inspiral of binary compact stars terminates by either 
the orbital instability (turning point of -Ebind) or the mass-shedding limit. In 
both cases, this defines the innermost stable circular orbit (ISCO). The end of 
inspiral phase strongly depends on EOS - for irrotational neutron star binaries 
a quasi-equilibrium sequence terminates by mass-shedding limit (circles at the 
end of each line in Fig 2) and for strange stars by orbital instability (shown as 
diamonds in Fig. 2). The differences in the evolution of binary strange stars 
and neutron stars stem from the fact that strange stars are principally bound 
by additional force than gravitation: the strong interaction between quarks. 
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They are self-bound objects having very high adiabatic index at the stellar 
surface (see (2q). Although the crust of a 1.35 M Q neutron star contains only a 
few percent of the stellar mass, this region is easily deformed under the action 
of the tidal forces resulting from the gravitational field produced by the com- 
panion star. The end of inspiral phase of binary stars strongly depends on the 
stiffness of matter in this region. The frequency of gravitational waves at the 
ISCO is one of potentially observable parameters by the gravitational wave 
detectors. In addition so called "break frequency", a characteristic frequ ency 



where the power emitted in gravitational waves decreases measurably [14j, [22 



could be an observable quantity. The frequencies of gravitational waves at the 
depature point from the 3PN aproximation and at the ISCO strongly depend 
on the EOS. For irrotational equal mass (of 1.35 M & at infinite separation) 
binaries the ISCO frequency is ~ 1100 — 1460 Hz for strange stars described 
by the MIT bag model and between 800 Hz and 1230 Hz for neutron stars de- 
scribed by nuclear EOS (see also 0,0]). The 3rd PN approximations for point 
masses derived by different authors are giving; ISCO at very high frequencies 
of gravitational waves > 2 kHz (0, 0, EH). 



It should be mentioned that the frequency at ISCO found in quasi-equilibrium 
approximation could differ from that obtained by the numerical integration 
of the full set of time dependent Einstein equations. The only comparison 
(according to the authors' knowledge) between these two approches for binary 
neutron stars can be found in Marronetti et al. (2004) who showed that for a 
polytropic EOS with 7 = 2 and M/R = 0.14 the frequency at the ISCO at 
quasi-equilibrium approximation is ~ 15% higher than the dynamical ISCO. 
Thus full relativistic hydrodynamical calculations are required in order to 
obtain detailed predictions for particular forms of the nuclear matter equations 
of state. 
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